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We examine in detail the changes in the production of primordial he- 
lium resulting from nonthermal neutrino momentum distributions produced 
by resonant transformation of electron-type neutrinos to steriles. These trans- 
formations, v e — > v s {y e — » z/ s ), amplify a positive (negative) lepton number 
asymmetry. We find that the resulting supression relative to a thermal dis- 
tribution of low energy v e reduces n — > p conversion to a greater extent than 
does the enhancement of n — > p from an identical suppression of D e . Thus, 
equal lepton-number asymmetries of opposite sign have unequal effects on the 
resulting helium yield in primordial nucleosynthesis. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35. +c; 95.30.-k 



1 



I. INTRODUCTION 



The role of neutrinos in primordial nucleosynthesis has been exploited to infer constraints 
on the number of leptons and light neutrinos in the Standard Model of particle physics [|l| . 
Complimentary to this, the width of the Zq was found to allow only three light weakly- 
coupled neutrinos 0. This experimental determination of the number of neutrinos, N u , has 
led some to describe primordial nucleosynthesis as determined by only one parameter: the 
baryon-to-photon ratio, rj. Along with the assumptions of homogeneity and isotropy, the 
theory of Standard Big Bang Nucleosynthesis (SBBN) predicts the abundance of the lightest 
nuclides within reasonable error through the single parameter rj. 

Now, however, indications of neutrino oscillations from the atmospheric muon neutrino 
deficit, the solar neutrino deficit and the LSND excess may prove that our understanding of 
neutrino physics is incomplete [[| . This uncertainty in neutrino physics threatens to destroy 
the aesthetics of a one parameter SBBN with the several parameters of the neutrino mass 
and mixing matrix. 

Neutrino mixing in the early universe can induce effects including matter-enhanced neu- 
trino flavor transformation, partial thermalization of extra degrees of freedom, and lepton 
number generation. Efforts are sometimes made to crowd the multiple effects of neutrino 
mixing into a single parameter, N^ G , the effective number of neutrinos. 

The kinetic equations describing the effects of neutrino mixing in the early universe were 
first discussed by Dolgov Q. Larger 5m 2 and sin 2 29 accelerate population of sterile states. 
Resonant conversion of active neutrinos to sterile neutrinos may also take place when the 
mass eigenvalue corresponding to u s is smaller than that for the u a , i.e., m 2 s —m 2 a = 5m 2 a < 0. 
Both resonant and non-resonant conversion has been used to place constraints on active- 
sterile neutrino properties 0-0]. 

Several formulations of neutrino dynamics in the early universe evolve the system at a 
monochromatic average neutrino energy with interactions occuring at an integrated rate. 
This approach is appropriate when the time scale of interactions between the neutrinos and 
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the thermal environment is much shorter than the dynamic time scale of neutrino flavor 
transformations. However, when tranformations and the Hubble expansion are more rapid 
than thermalization, neutrino energy distributions can be distorted from a Fermi-Dirac 
form. Furthermore, such spectral distortions can survive until nucleosynthesis ||. The role 
of neutrino energy distributions in primordial nucleosynthesis has been discussed since the 
earliest work on the subject |9,10|]. Recent work |Tl~l , |l2ll has emphasized the importance of 



neutrino energy distributions on the weak nucleon interconversion rates (n ^± p). Kirilova 
and Chizhov [|llj included the distortion of neutrino spectra from nonresonant transforma- 



tions in their nucleosynthesis calculation. A time-evolving distortion was also included in 
a calculation of helium production after a neutrino-mixing-generated lepton asymmetry in 
electron- type neutrinos from v e ^± v s or v e ^± v s (active-sterile) mixing |12| . 



In this paper we explore in detail the effects on primordial nucleosynthesis of the non- 



thermal neutrino spectra produced by the resonant transformations discussed in The 
resonant transformation of electron-type neutrinos to steriles can leave a near absence of 
neutrinos in the low energy end of the active neutrino energy spectrum. We find that this 
suppression of low energy neutrinos or anti-neutrinos in the otherwise Fermi-Dirac spec- 
trum alters the neutron-to-proton weak interconversion reactions. These alterations stem 
from initial and final phase space effects. 

Specifically, we find that the absence of low energy electron anti-neutrinos lifts the Fermi 
blocking of final state neutrinos in neutron decay. In addition, the absence of low energy 
electron neutrinos capturing on neutrons significantly alters the rate of this interaction. 
Both cases lead to an alteration of the neutron-to-proton ratio at weak freeze out, which is 
directly related to the produced mass fraction of primordial helium. 

II. THE DYNAMICS OF NEUTRINO TRANSFORMATION 

The general equations describing the evolution of M flavors of neutrinos in a high density 
environment such as the early universe are derived from the evolution equation for the density 
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operator 



i% = \HM. (1) 



Here p can be the density operator for the entire system of M neutrinos or any smaller 
neutrino system that evolves rather independently of the other neutrino mixings. Several 
authors explicitly follow the evolution of neutrinos through the density matrix. 

The evolution of the density operator for a two state active-sterile system has been 
popularized in a vector formalism where the density operator is 

p = ip (l + p.a) p=ip (l + P-a) (2) 

where p, P and Pq refer to the anti-neutrino states. For an active-sterile two state system, 
the Hamiltonian of eqn. (|l|) includes forward scattering off of the background, and elastic 
and inelastic scattering of neutrinos from all leptons (see, e.g., Ref. [13]). For this system, 
substituting @ into the evolution equation resulting from Eqn. [I] gives 

ip = <V> X P + (1 - P.) (j t In P ) z - (d e + D 1 + j t In P ) (P,x + P y y) (3) 



dt 

d 



P = ( T ( u eVe -> aa))(\ a n a n a - n u n Ve ) (4) 



Qt=eM a M T 



where an average over the momentum ensemble of neutrinos is taken for the vector potential 
V and the interaction rates V. Here D 1 and D E are the damping coefficients due to inelastic 
and elastic neutrino scattering off of and themselves (see 0). Also, \ v = \ and A e = 1. 
This set of equations (HfD can be followed for specific neutrino momenta. An initial thermal 
neutrino spectrum can be discretized into M energy bins. Along with the number density 
evolution equations of neutrinos that do not undergo transformation, these equations form 
8J\f equations describing the evolution of the system. 

Here we assume that the sterile neutrino sea is unpopulated at a temperature T ~ 
100 MeV (see Ref. ||14j| ). Active neutrinos can undergo resonant transformation to sterile 
neutrinos as the density of e ± decreases with the expansion. As in the sun, the resonance 
point is energy dependent. For example, the resonance, for 5m 2 sa < 0, occurs at 



E\ \6m?/eV 2 



(5) 



. ■>- / res 16(T/Mev) 4 L„ Q ' 

where E/T is the neutrino energy scaled by the ambient temperature. This resonance 
adiabatically converts v a (i/ a ) to v s {j? s ) as long as the timescale of the resonance is much 
longer than the neutrino oscillation period at resonance, and the resonance position must 
move slowly through the neutrino energy spectrum (see Ref. fl2fl). 

As the universe expands and the ambient temperature decreases, the resonance moves 
from lower to higher neutrino energies. The energy-dependence of the resonance becomes 
important at temperatures below T^ ec « 3MeV ||. Below this temperature, the neutrinos 
decouple from the background. The weak interaction rates of elastic and inelastic scat- 
tering of neutrinos above Td ec are much faster than the expansion time scale. In this case, 
the neutrino energy distributions can reshuffle into a thermal Fermi-Dirac spectrum before 
the onset of nucleosynthesis. However, if the transformation occurs after decoupling, the 
expansion rate will be too rapid, and therefore the neutrino spectra enter the epoch of nu- 
cleosynthesis in a non-thermal state. This will be the case, usually, when the mass-squared 
difference between mass eigenvalues corresponding to v e and v s is |<5mgJ < leV 2 . 



cLS cL 



Direct v e — > v s {p e — > v s ) transformation was explored by Foot and Volkas [15 
possible means to allow for a sterile neutrino within SBBN with minimal modification to 
the theory, and also to reduce the difference between the 4 He abundance predicted by the 
observed primordial deuterium abundance jn[ and the primordial 4 He abundance inferred 
by Olive et al. Jl6|]. The change in primordial helium production in this v e — > v s {p e — > i> s ) 



scenario, including alterations of neutrino spectra, was described in [0. The resonant 
transformation of u e — > u s (9 e — > v s ) produces a lepton number asymmetry in electron- 
type neutrinos, L Ve . Here we define the net lepton number residing in neutrino species 
u a to be L Va = {n Va —n Pa )/n^, where n Va , n Va , and n 7 = (2((3)/tt 2 )T 3 are the proper 
number densities of u a , u a , and photons, respectively. The sign of this asymmetry depends 
sensitively on initial thermal conditions and can be different in causally disconnected regions 
Hl8|Jl9|| . In the following sections, we address how non-thermal spectra affect the production 



of primordial helium. 



III. WEAK NUCLEON INTERCONVERSION RATES AND PRIMORDIAL 

HELIUM 

Since nearly all of the neutrons present during primordial nucleosynthesis go into 4 He 
nuclei, the production of primordial helium in the early universe is dominated by the abun- 
dance of neutrons relative to protons. The weak nucleon interconversion reactions 

n + v e p + e~ n + e + ^ p + u e n ^ p + e~ + v e . 

are in steady state equilibrium at temperatures T > 0.7 MeV. The neutron-to-proton ratio 
(n/p) freezes out of equilibrium when the weak nucleon interconversion rates fall below 
the expansion rate (see Fig. [I]). The "frozen" n/p ratio slowly evolves as free neutrons 
decay until the epoch (T ~ 0.1 MeV) when almost all neutrons are incorporated into alpha 
particles ("nucleosynthesis"). In steady state equilibrium, n/p can be approximated by the 
ratio of the nucleon interconversion rates, 

1 

P A n ^p 

The evolution of n/p is followed numerically in detail in Kawano's J20| update of Wagoner's 



code ||10|| . The suppression of low energy electron- type neutrinos affects the rates of the six 
weak nucleon reactions differently. 

If sterile neutrinos are partially or completely thermalized, they will increase the expan- 
sion rate since they will contribute to the energy density of the primordial plasma. The 
expansion rate, H, is related to the total energy density, p to t and the cosmological constant, 
A as [0,0] 



3 V 3 . 

In general, an increase to the energy density and expansion rate will increase n/p at freeze- 
out, and thus increase the primordial 4 He yield, Y„. The contribution to the expansion rate 



is not significant for v e — *> v s [p e — > z/ s ) conversion post-decoupling, since the energy density 
in electron-type neutrinos is just transferred to steriles without further thermal creation of 
neutrinos (since they have decoupled from the plasma before the onset of transformation). 

For the case where L Ue < 0, the reactions influenced by the v e distortion are the forward 
and reverse reactions of n + v e ^ p + e~. Because of the neutron/proton mass difference, 
very low energy neutrinos readily participate in n + v e — > p + e~. Here, the neutrino 
energy is E v = E e — Q > 0, where Q = m n — m p is the nucleon mass difference, and 
E e is the electron energy. The suppression of low energy neutrinos impacts this reaction by 
significantly limiting the phase space of initial states of the interaction. The neutrino energy 
distribution is accounted for in the n + v e — > p + e~ rate integral in the following manner: 

Xn^ pe = Aj v e Ey u dp u [e E »/ kT » + l]- 1 ^ + e- E °l kT ]-\ (6) 

The rate of the reaction p + e~ — > n + v e is given by 

A P ™ = A J Ey e dp e [e E ^ kT + l^l + e - E ^]-\ (7) 

In these equations we have used Weinberg's notation [21]. The deficit in low energy neutrinos 
renders the integrand of (|5]) to be nearly zero for energy ranges where v e have transformed 
to steriles (see Fig. ||]). The reaction p + e~ — > n + v e is altered by a lifting of Fermi- 
blocking [1 + g-^A^]- 1 a t low energies because of the reduction of low energy neutrino 
numbers. However, the rate \ pe ^ nu is changed less significantly than X nu ^ pe and thus affects 
n/p less. This can be seen in Fig. ^ (a), where the change in rate with respect to the overall 
n — > p{p — > n) rate is shown for rate (2) \ nu _, pe and (6) X pe ^ nu . 

For L Ue > 0, the reactions affected by i> e distortion are n+e + ^ p+v e and n ^ p+e~ +i> e . 
In the forward and reverse reactions n + e + ^ p + u e , the escaping or incident v e must have 
an energy E Pe = Q + E e > 1.804 MeV. Even for the case we examined with the greatest 
spectral distortion (|5m^ e | = 1 eV 2 ), the spectral cutoff never extends above E Pe = 1.35 MeV. 
Therefore, n + e + ^ p + v e is not altered by this low energy distortion (see Fig. ^ (b)). The 
interaction p + e~ + v e —>■ n makes a very small contribution to X p ^ n , and modifications of 



its rate are not important to n/p at freeze-out. The only significant reaction rate affected 
by the v e distortion is neutron decay: n — > p + e~ + v e (Fig. ^). The reaction is limited in 
SBBN by Fermi-blocking of low-energy v e products. This blocking is lifted when low-energy 
v e transform to v s . The integrand of the neutron decay rate 

K-> P ev = Aj v e ElE 2 e dp v [l + c -&/«^]-i[i + c -&/wv]-i (8) 

can be seen in Fig. [|for the most dramatic case of \5m 2 \ = 1 eV 2 . The Fermi-blocking term, 
[1 + e - E "/ kT »]~ 1 ^ i s unity for low energy neutrinos that have transformed to steriles. 

The change in \ nv ^ pe (for L Ve < 0) is significantly larger than that for X n ^ pev (for 
L Ve > 0). The change in helium mass fraction (Y p ) produced by BBN for L Ue < is 



SY p ~ 0.022, while a L Ue > correspondingly causes smaller a change of 5Y P ~ —0.0021 |T2 

In summary, if electron-type neutrino spectra are altered through matter-enhanced trans- 
formation and survive until nucleosynthesis, they can modify n ^ p interconversion rates 
through the availability of initial and final neutrino energy states. This kind of scenario is 
realized in u e — > v s {p e — > u s ) matter-enhanced transformations with \Sm 2 \ < leV 2 . Non- 
thermal v e j v e spectra will arise in any case where flavor transformations occur near or below 
T dec between a more-populated electron-type neutrino and some other less-populated neu- 
trino flavor, and will produce the same kinds of limitations on reaction phase space as the 
resonance point goes from lower to higher neutrino energies. The effects of non-thermal 
electron-type neutrino distributions will then need to be included. 
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FIG. 1. The evolution of the n f± p rates are shown. 1) n — > p + e~ + P e , 2) n + v e — > p + e~ , 
3) n + e + — > p + v e 4) p + P e — > n + e + , 5) p + e~ + v e — > n, 6) p + e~ — > n + v e 7) The Hubble 
expansion rate, when the rates fall below the expansion rate, n/p "freezes" out. 



a b 




4 6 8 10 4 6 8 10 



E„/rri E /rn 

e' e e' e 

FIG. 2. Plotted are the contours of the amplitude of the rate integrand for n + v e — > p + e _ , 
with increasing amplitude towards the upper right: (a) shows the blocking of the integrand for 
lower energies due to the suppression of low energy electron neutrinos; (b) shows the standard 
BBN integrand. 
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FIG. 3. Plotted are the contours of the amplitude of the rate integrand for n — > p + e~ + z/ e , 
with increasing amplitude towards the center, (a) shows the lack of Fermi blocking by low energy 
v e at higher temperatures, (b) shows the standard BBN integrand. 
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FIG. 4. Plotted are a) the change in rate 2. n + v e 
n — > p + e~ + u e , 3. n + e + — > p + P e 4. p + P e — > n + e H 
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p + e~, 6. p + e _ — > n + v e \ and b) 1. 
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